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ABSTRACT 

We study the role of the unknown microphysical properties of carbonaceous dust par- 
ticles in determining the amount and size distribution of carbonaceous dust condensed 
in type II supernova explosions. We parametrize the microphysical properties in terms 
of the shape factor of the grain and the sticking coefficient of gas-phase carbon atoms 
onto the grain surfaces. We find that the amount of dust formed is fairly independent 
of these properties, within the parameter range considered, though limited by the 
available amount of carbon atoms not locked in CO molecules. However, we find that 
the condensation times and size distributions of dust grains depend sensitively on the 
microphysical parameters, with the mass distributions being weighted toward larger 
effective radii for conditions considering grains with higher sticking coefficients and/or 
more aspherical shapes. We discuss that this leads to important consequences on the 
predicted extinction law of SN dust and on the survival rate of the formed grains 
as they pass through the reverse shock of the SN. We conclude that a more detailed 
understanding of the dust formation process and of the microphysical properties of 
each dust species needs to be achieved before robust prediction on the SN dust yields 
can be performed. 
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1 INTRODUCTION 

Interstellar dust, once considered to be little more than 
a nuisance to astronomical observations, is one of the 
most interesting area s of astrophysical research today 
(|Li fc Greenbera 120031 ). One aspect of particular interest 
is the origin of dust at high redshift [z > 5). Possible 
sources of interstellar dust that have been considered 
include outflows from a sympt ot ic giant branch (AGB) stars 
dMorgan fc Edmunds! 120031 ; IZhukovska. Gail, fc Trielofl 
20081 : IValiante et alj|2009h. Wolf-Rayet systems llCherchnefl 
201(j ), quasars (|Elvis. Marengo, fc Karovska 2002) , and 
supernova explosions (jKozasa, Hasegawa. fc Nomotol 
19891. Il99j; ITodini fc Ferraral l200ll; INozawa et all 120031; 
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possible sources at such high redshift remains unclear. 
Supernovae (SNe) are considered by some to be likely 
contributors of much of the dust in the early Universe 
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have not been observed, making determinations of their 
dust contributions difficult. Instead we need to consider 
dust yield predictions of more recent SNe, which have 
been observed, and then extrapolate the dust yields to 
earlier SNe. The the oretical predictions of dust yie l ds for 
recent supernovae iKozasa. Hasegawa. fc Nomotol Il99ll : 
ITodini fc Ferraral l200ll : INozawa et al.ll200Sl) . however, are 
too large compared to observations, with discrepancies 
that ca n be as high as 3 to 4 orders of magnitude in some 
cases llHovle fc Wickramasinghel Il97d; iLucv et all Il989l: 
" 1 19931: lElmhamdi et all 120031 ; iMeikle et all 
Kotak et al.ll2009i ). Observations of young supernova 
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remnants (SNR), however, have confirmed dust masses one 
to two orders of magnitude greater than SNe. Some of these 
SNR have observe d dust masses of 0.02-0.05 4 Ma by Spitzer 
llRho et alJ I2008D . 0.06 M Q by A KARI llsibthorpe et ail 
|2010() , and 0.075 M Q by Hershel (|Barlow et alj|2010h . in 
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Cassiopeia A, and 0.04- 0.1 Mp by Spitzer in the pulsar 
wind nebula G54. 1+ 0.3 (|Temim et alj|20ich . Additionally, 
iNozawa et al . (2010) have demonstrated that the observed 
spectral energy distribution of Cassiopeia A can be well 
reproduced by the calculations of dust formation in Type 
lib SNR and the mass in the SNR is 0.07 M s . 

Reconciling the dust yield prediction for local SNe can 
be done via one of two channels: improving our understand- 
ing of the dust formation process, leading to a substan- 
tially decreased prediction, or revising the observational con- 
straints to account for a higher dust yield than so far im- 
plied. It is possible that some amount of dust has avoided ob- 
servation. Dust particles absorb light and re-emit the energy 
at infrared wavelengths. However, cold dust at temperatures 
of a few tens of Kelvin could escape detection at mid-infrared 
wavelengths. Additionally, areas where the dust is optically 
thick could obscure some amount of dust, again allowing 
some dust to not be detected. Dust clumpin g may also affect 
the e stimates of dust mass from absorption (|Sugerman et alj 
2006), since it is usually assumed that the surface filling 
factor of dust is close to unity, while substantial clump- 
ing could be present due to the intrinsic inhomogeneity of 
the eje cta and Ray leigh- Taylor instabilities in the expandin g 
eject a l|Woodedll9lmlDouvion. Lagage. fc Cesarskvlll999l ). 

On the theoretical side, all the estimates for dust pro- 
duction in SNe ar e based on the so-called classical theory 
of nucleation (e.g.. [Becker fe Dorind |l935l )). It has been ar- 
gued, however, that the use of classical nucle ation theory in 
astrophysical enviro nments is questionable (jPonn fc Nuthl 
1 19851 ; lLazzatilliooil '). I n addition, almost all SN dust nu- 
cleation models thus far have considered the formation of 
spherical grains, and assumed any atoms/molecules that 
contact the grain will adhere to the grain; conditions that 
reflect maximally efficient nucleation. It is therefore not en- 
tirely surprising that the theoretical estimates of SN dust 
yields, based on upper limit of nucleation efficiency, are in 
excess of those from observ ations. It should be noted that 
iBianchi fc Schneider! (|2007l ) have considered less than max- 
imally efficient nucleation by assuming the probability of 
atoms/molecules adhering to the grain is less than unity, 
resulting in smaller dust yields. 

Since dust formation is a highly non-linear phe- 
nomenon, understanding the effects of different nucleation 
rates on the final dust yields is difficult. To check the ef- 
fects of different nucleation conditions we have performed a 
parametric study, in which we consider carbonaceous dust 
production in a SN explosion by varying the shape of the 
forming grains as well as the sticking coefficient, i.e., the 
probability that an incoming monomer will stick to the 
grain rather than bounce off and remain in the gas phase. 
Our study is phenomenological and aims at understanding 
which conclusions of previous nucleation studies are robust 
to changing the parameters, and which may need to be 
investigated more thoroughly. A self-consistent nucleation 
model in astrophysical c onditions will be achieved by involv- 
ing a kinetic approach (|Donn fe Nuthl 1 19851 ; lLazzatil 120081 ; 
iKeith fc La zzati 20 11]) and detailed chem istry of precursor 
molecules j|Cherchneff fc Dwek|[2009l . l2010l ). Such a detailed 
approach, however, is still under development and is not yet 
applicable to large scale simulations like the one that we use 
here , and that have been us e d in previous investigations of 
e.g jTodini fc Ferraral (|200ll ); INozawa et all (|2OOa . l2O10h . 



This paper is organized as follows: in Sect. 2 we detail 
the nucleation theory that we adopted; in Sect. 3 we describe 
the numerical code used for the computations; and in Sect. 4 
we describe our results. In Sect. 5 we finally discuss the 
implication and limitations of our results and lay out future 
perspectives for SN dust studies. 



2 NUCLEATION 

Nucleation is the first step of a first-order phase transition. 
In the case we consider here, the phase transition is from a 
gas of carbon atoms to solid clusters of amorphous carbon 
dust. There is a phase equilibrium pressure where both the 
gas and solid phases are stable within a given volume. The 
phase equilibrium pressure depends on the temperature of 
the materials in the gas phase within the volume. Nucleation 
of clusters of atoms/molecules of the new phase is favoured 
when the gas phase is supersaturated (i.e., the pressure of 
the gas is higher than the phase equilibrium pressure). The 
higher the supersaturation (i.e., the ratio of the pressure 
to the phase equilibrium pressure), the smaller the size of 
the stable clusters that are able to form. The size of the 
smallest stable cluster able to form at a given temperature 
and density is called the critical cluster size. Clusters that 
are smaller than the critical size will tend to evaporate, while 
larger clusters will tend to continue to grow. The goal of any 
nucleation theory is to calculate how many critical clusters 
form per unit volume per unit time. 

The classical theory of nucleation considers nucle- 
ation as a thermodynamical process in quasi-equilibrium 
( Bec ker fc Dorind[l935l : iFeder et al ] |l966l : iKashchievllioOol ) . 
Besides the supersaturation level, the physical properties 
of the nucleating material affect the size of a critical clus- 
ter. The classical nucleation theory assumes that all clus- 
ters share the same properties, such as surface tension and 
shape, independent of their size. The theory also assumes 
that those properties are equal to those of a macroscopic 
sample. Moreover, the clusters are assumed to have a uni- 
form equilibrium temperature that is equal to the temper- 
ature of the surrounding gas. These basic assumptions are 
problematic because macroscopic thermodynamic properties 
are not expected to be applicable to clusters of only a few 
atoms. 

A different approach is provided by the kinetic theory 
of nucleation, which is applicable to very small cluster sizes. 
The kinetic theory relies on calculating the attac hment and 
detac hment frequencies of monomers to a cluster (Kashchicv 
2000). Furthermore, in contrast to the classical nucleation 
theory, the kinetic theory follows the formation of clusters 
smaller than the critical cluster size. In the framework of the 
kinetic theory the critical cluster is the cluster whose attach- 
ment and detachment frequencies are equal, and thus it is 
stable. The downside of the kinetic theory is that attach- 
ment and detachment frequencies for all cluster sizes need 
to be calculated in order to determine the overall nucleation 
rate. The main aim of this paper is to reveal the dependence 
of dust formation processes on the microphysical properties 
of grains. In order to achieve this, we adopt the simpler clas- 
sical nucleation theory, rather than using the kinetic theory 
that demands a more complicated treatment. 

When the supersaturation level is 5* > 1, nucleation can 
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take place because the free e nergy of the new phase is lower 
than that of the old phase (|Kashchievll200(]h . The change 
in the free energy is due to the work necessary to form the 
critical size clusters. Energy is released in the formation of 
the volume of the cluster, but ener gy is required in order 
to form the surface of the cluster jKashchievll2000l ). The 
nucleation rate (Eq. (HI) is given basically by two factors: the 
number density of critical clusters and a kinetic factor that 
describes the rate at which clusters become large enough to 
be stable (i.e., critical size or larger). 

The general eq uation of stationary nucleation is given 
by (jKashchiev! (|2OO0h . their equation 13.39): 



Js = Aexp(-W*/kT), 



(1) 



where A is the kinetic factor, W* is the work needed to 
form the critical size cluster from the gaseous state, k is 
the Boltzmann constant, and T is the gas temperature. The 
kinetic factor, A, is given by A = A' exp (A/i/fcT), and A' 
is: 



A = 7 



187r 2 mo 



\ kT ) 



Co, 



where 7* is the size dependent sticking coefficient, c is 
the shape factor of the cluster, a is the surface tension, 
mo and vq are, respectively, the mass and volume of the 
monomer of nucleating material, p e is the phase-equilibrium 
pressure, and Co is the concentration of sites where nucle- 
ation can occur (equations 13.41 and 13.44, respectively in 
lKashchievll|2000t )'). The concentration of gaseous monomers, 
Ci is related to the concentration of nucleation sites by 
Ci = Co exp ( — Wi/kT), where Wi is the work needed to 
form a cluster co nsisting of one monomer (equation 7.5 in 
lKashchievi(|2000l )'). We consider the monomer in the gaseous 
state to be indistinguishable from the monomer in the con- 
densed state, so that Wi = and Ci = Co. Using the su- 
persaturation ratio, 5* = Ci/Ci, e = p/Pe and A/j, — kTlnS, 
where Ci, e = p e /kT is the gaseous monomer concentration 
at the phase-equilibrium pressure and p is the partial pres- 
sure of gaseous monomers, we obtain 



A = j' 



( c 3 v 2 a 
\ 187r 2 mo 



While the sticking coefficient may depend on the size of the 
cluster, we consider it to be constant, so that 7* = 7 = con- 
stant. Finally, using th e work to fo r m a c ritical sized cluster 
W* = 4c 3 'via 3 1 '27 A/x 2 (|Kashchievl (|2000T ). equation 4.8), we 
find our stationary nucleation rate equation to be: 



J s = 7 



3 2 
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(2) 



After nucleation the clusters grow through impingement 
of monomers upon the cluster. To find how much the clusters 
grow we begin by finding the volume of the newly nucleated 
cluster. The clusters nucleate with some critical number of 
monomers, n* (equation 4.7 in iKashchievI (|2000h ): 



8c 3 v 2 a 3 



(3) 



27A M 3 ' 

The volume of the critical cluster is then, v* = n*vo. The 
change in the cluster's volume over time depends on the 
sticking coefficient 7, the surface area of the cluster S, the 



concentration of monomers Ci , the volume of the monomer 
vo, and the average relative speed of the monomers with 
respect to the cluster. In this paper we consider clusters 
that can nucleate with aspherical shapes through the use 
of the shape factor c. The shape factor is a dimensionless 
quantity that relates the su rface area S of a n object to its 
volume V by: c = E/V 2/3 jKashchiev! l200oh . Thus we find 
the change in volume over time to be: 



dt \ 2-KTTlo J 



(4) 



In the same manner as iNozawa et ail (12003T ) we compute 
the depletion of the available nucleation material through 
mass conservation: 



Ci(t) 
Ci(t) 



1-Vi 



J{t') V(t,t') 
Ci(t') 



I'll 



dt', 



(5) 



where Ci is the nominal concentration of monomers - the 
concentration expected should nucleation not occur, which 
in an expanding shell of volume V s heii can be found using 

Cl(t = t n )Vshell(t = t n ) = C\(t = to)V s hell(t = to), 

so that the total number of gas-phase atoms remains con- 
stant, and V(t, t') is the volume of a cluster formed at time 
t' and measured at time t. Rather than computing the in- 
tegral on the right-hand side, we instead calculate Ci(t) as 
described in step (iii) in the next section. From Eqs. ((2)1 
and Q, we see that J s and dV/dt are simply proportional 
to the sticking coefficient 7. Thus, we expect that reduced 
sticking coefficients will suppress both nucleation and grain 
growth. The dependence of J s on the shape factor c is more 
complicated, since it appears in both the kinetic factor and 
the exponential term. The shape factor in the exponential 
term, however, will dominate the nucleation rate equation 
and we expect that increased shape factors will suppress nu- 
cleation. On the other hand, grain growth (dV/dt) is simply 
proportional to c, and thus an increase in the shape factor 
will increase the cluster growth rate. 



3 SIMULATIONS 

We concentrate on the formation of carbonaceous grains 
(clusters) from carbon atoms in the expanding material of 
a core-collapse supernova. Our simulations are based on the 
hydrodynamic results and elemental composition for the 
unmixed ejecta of a core-collapse supernova of a 20 Mq 
progenitor star with metallicity Z = and an explosion 
energy of 10 51 ergs by lUmeda fc Nomotol l|2002l 'l (see also 
iNomoto et all |2006l )). 

Table [T] describes the data necessary for the calculation 
of carbon grain formation. To compute the supersaturation 
of the expanding gas, we find the phase equilibrium pressure 
by: p e — poe~ A / T+B , where pa is the standard pressure, T 
is the temperature, and the values of A and B, list e d her e 
in Tabled] are taken from Table 2 of INozawa et al.l (|2003h . 
As the material ejected by the core-collapse supernova ex- 
pands it also cools. Figures Q] and [2] show the evolution of 
the density and temperature, respectively, of the expanding 
and cooling material for two of the enclosed mass subshells 
that we refer to in the rest of this work. 

We study four different values of the unknown sticking 
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Figure 1. Density evolution for 4.96 (black) and 6.0 (red) Mg 
enclosed mass subshells up to 1000 days after the SN explosion. 
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Figure 2. Temperature evolution for 4.96 (black) and 6.0 (red) 
Mq enclosed mass subshells up to 1000 days after the SN explo- 
sion. 



coefficient, 7 = 1.0, 0.1, 0.01, and 0.001, neglecting any de- 
pendence of 7 on the temperature of the gas and the size 
of the cluster. For each value of the sticking coefficient, we 
study six different values for the shape factor, c = (367r) 1/ ' 3 , 
5.4, 6.0, 7.0, 9.0, and 12.0, corresponding to shapes rang- 
ing from a sphere to a flattened cylinder similar to a coin. 
We therefore performed 24 simulations in total. A sticking 
coefficient of 7 = 1.0, and a shape factor for a sphere of 
c = (367I") 1 / 3 are the usual parameters used in previous nu- 

cleation w orks. 

As in iNozawa et alJ <|2003l ), we assume that the stable 
formation of CO molecules occurs prior to grain nucleation 
and that carbon grains will form only where the initial num- 

Table 1. Carbon properties 



A/10 4 
(K) 



(erg-cm 2 ) 



VQ TOO 

(10- 24 cm- 3 ) (10- 23 



Figure 3. Number fraction of helium (green dash-dot), carbon 
(black dashed), oxygen (red dotted), and silicon (magenta dash- 
dot) atoms for enclosed masses from 2.45 to 8.7 Mq. The solid 
(blue) line is the carbon number fraction after the formation of 
CO molecules. 



ber fraction of carbon is higher than that of oxygen. Under 
this assumption, the number fraction of free carbon atoms 
available for dust formation is obtained simply from the 
initial number fraction of carbon minus the number frac- 
tion oxygen. We divide the expanding gases into a series of 
enclosed mass shells beginning at ~ 4.93 Mq and ending 
at ~ 6.21 Mq. In this range of enclosed masses the num- 
ber fraction of carbon atoms, after the formation of CO, 
is highest; ranging between 2 x 10 _1 and 8 x 10" 9 . Fi euro 
[3] shows the number fraction of carbon and oxygen atoms 
in the expanding eject a from the hydrodynamic results of 
lUmeda fc Nomotol (|2002h . The solid (blue) line indicates the 
number fraction of carbon left over after the formation of CO 
molecules. It should be noted for completeness that grain 
nucleation can occur at enclosed masses larger than we con- 
sider here, but is extremely inefficient due to low carbon 
number fractions, ~ 10~ 14 and below. 

For each mass subshell, our code starts by following the 
evolution of the density and temperature (see Figures [T] and 
[2] respectively, for examples) of the gas until the condition 
of supersaturation is satisfied. From that point on, at each 
time step the code performs three operations. 

(i) First, the code computes the number of critical clus- 
ters that are formed given the saturation, temperature, and 
partial pressure of the carbon atoms (according to Eq. @). 

(ii) Second, the code grows any pre-existing grain formed 
at earlier times according Eq. Q. 

(iii) Finally, the code subtracts from the carbon in the 
gas phase the amount of carbon that has been locked in 
the solid phase by the processes in steps (i) and (ii). The 
concentration of gas-phase carbon is evolved according to: 



Ci(*n) 



where AV a < 



Ci(t n -i) — 

Cl(tn-l) 

Ci(t n ) 



•ins, 11 
VoV s hell,n 

(&(t„-i)- 



Ci (*«)), 



8.64726 19.0422 



1400 



8.805 



1.995 



n is the total change in volume of grains from 
the previous time step to the current time step. Then, the 
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Figure 4. Nuclcation rates for spherical grains ( c = (367T) 1 / 3 ), 
for the four considered sticking coefficients of 7 = 1.0, 0.1, 0.01, 
and 0.001. Filled circles indicate the time and rate of maximum 
nucleation. Rates are calculated at an enclosed mass coordinate 
4.96M Q . 



first term on the right-hand side is the concentration of gas- 
phase carbon from the previous time step, the second term 
is the change in the concentration of solid-phase carbon (it 
is a summation of all grain changes, including formation of 
new critical clusters and the growth of existing grains), and 
the final term accounts for the decrease in concentration of 
gas-phase carbon due to the expansion of the shell. 

This process is repeated for each mass shell until the con- 
centration of gas-phase carbon is reduced to 1 per cent of 
its original value. The dust distributions from all shells are 
then summed together to produce the final dust yield of each 
particular set (7,0). 
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Figure 5. Maximum nucleation rates for spherical carbon grains 
at enclosed masses < 6.2 Mq for four sticking coefficients. 
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4 RESULTS 

4.1 Spherical Grains 

After the supernova explosion, the hot gases expand and 
cool. The cooling allows the gases to reach supersaturation 
conditions. Once the supersaturation level becomes greater 
than unity, nucleation can occur. The supersaturation level 
continues to rise and the nucleation rate increases over time 
until depletion of available material becomes significant and 
the supersaturation level begins to drop, after which the 
nucleation rate falls off quickly. Eventually the gas is no 
longer supersaturated and nucleation ceases. However, grain 
growth is still possible. 

Figure [4] shows the nucleation rate for spherically 
shaped grains at enclosed mass coordinate 4.96 M© for the 
four sticking coefficients. We chose this particular enclosed 
mass shell because it contains the highest abundance of car- 
bon atoms, after CO formation, of all our mass shells. The 
solid curve corresponds to c = (367r) 1//3 and 7 = 1.0, the pa- 
rameters generally used for nucleation studies. The dashed, 
dash-dot, and dotted curves correspond to 7 = 0.1, 0.01, 
and 0.001, respective l y. To be consistent with the works of 
IKozasa &c Hasesawal (I1987I - ) ; iKozasa. Hasegawa. fc Nomotol 



1 19891. Il99ll ); iNozawa et all l|2005t ). we consider the time at 



Figure 6. Nuclcation rates for spherical carbon grains at an en- 
closed coordinate of 6.00 Mq for four sticking coefficients. Con- 
densation times are indicated by filled circles. 



which the nucleation rate is at its maximum to be the con- 
densation time of the grains. We show these condensation 
times as filled circles at the maxima of the nucleation rates 
in the figure. 

At early times, the reduced sticking coefficient makes 
the formation of critically sized grains more difficult and 
results in a suppressed nucleation rate. In the absence of 
strong nucleation, carbon atoms are not depleted from the 
gas and the saturation continues to increase. The reduced 
sticking coefficient thus causes the time at which nucleation 
is at its maximum to be delayed, and the nucleation to take 
place at higher saturation levels. As a consequence, a larger 
number of critical clusters can form with much smaller size 
(Eq.EJ). 

Figure[S]shows the maximum nucleation rates for spher- 
ical grains as a function of enclosed mass for our four sticking 
coefficients. For enclosed masses up to M r ~ 5.87 Mq, the 
nucleation rate maxima for reduced sticking coefficients ex- 
hibit similar behaviour as in Figure [4] At greater enclosed 
masses, however, the behaviour is inverted and a reduced 
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Figure 7. Condensation times of spherical carbon grain as a 
function of enclosed mass for four sticking coefficients. 



sticking coefficient results in depressed maximum nucleation 
rates. This difference in behaviour is due to the reduction 
in the number fraction of carbon available, from ~ 10 -2 ' 5 
to ~ 10 -5 (see Figure [3]), which corresponds to a reduction 
in the concentration of carbon monomers. Since the nucle- 
ation rate (Eq. p]l) is proportional to the concentration of 
monomers squared, the drop in the carbon concentration 
consequently drops the nucleation rate. Thus, the nucleation 
rates do not peak as strongly (Figure |6| , drawing out nu- 
cleation to later times for all sticking coefficients, so that 
a catastrophic reduction in the available material does not 
occur. Thus, the nucleation rates for reduced sticking coeffi- 
cients remain depressed throughout the simulation time for 
higher enclosed mass shells. 

Figure[7]shows the condensation times corresponding to 
the maximum nucleation rates shown in Figure [5] The con- 
densation times of the solid curve (c = (367r) ly/3 , 7 = 1.0) 
are in good agreement with the condensati on times for car- 
bon grains reported in lNozawa et alJ (|2003h . Here the lower 
sticking coefficients result in delayed condensation times up 
to M r ~ 5.87 Mq, outside which the condensation time is 
only slightly delayed, or no longer delayed, when compared 
to the case of 7 = 1.0. The more noticeably delayed conden- 
sation times, as well as the similarity of condensation times 
for all our sticking coefficients, at higher enclosed mass shells 
are also due to the drawn out nucleation process already dis- 
cussed above. 

Reducing the sticking coefficient makes both nucleation 
and grain growth more difficult. A reduced sticking coeffi- 
cient results in a larger number of smaller grains, fewer larger 
grains, and a smaller maximum grain radius. In Figures [8] 
and [9] respectively, we show the size and mass distributions 
for spherical grains. For the case of 7 = 1.0, the maximum 
grain size achieved is between 2 and 3 /jm. When the stick- 
ing coefficient is reduced to 0.001, the maximum grain size 
is decreased to less than 0.01 fim. 

Since the reduced sticking coefficient causes larger num- 
bers of small grains to form, the small grains contain more 
relative mass than the larger grains, as can be seein in Fig- 
ure [9] However, the total mass of the grains is relatively 
robust. The total masses of dust grains for the spherical 



Figure 8. Size distribution of carbon grains for 7 = 1.0, 
0.1, 0.01, 0.001, with c = (367T) 1 / 3 . For reference, the solid 
(red) line represents the power-law distribution with the form 
of N r oc r ~ 3,5 , which has been suggested as th at of interstellar 
grains (e.g., Ma this. Rumpl. fc Nordsieckl ljl977h l. 
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Figure 9. Mass distribution of carbon grains for 7 = 1.0, 0.1, 
0.01, 0.001, with c = (367T) 1 / 3 . 

case are shown in Table [3] along with the aspherical cases 
which are discussed in the next section. Even though the 
onset of nucleation is delayed due to the reduced sticking 
coefficient at enclosed masses less than 5.87 Mq, where the 
majority of available carbon is contained, the subsequent 
grain growth consumes almost all of the carbon atoms for 
7 > 0.001. Thus, the total mass of carbon dust is principally 
determined by the mass of pre-existing carbon atoms. 



4.2 Non-spherical Grains 

We also calculated nucleation rates for aspherical grains. We 
choose a range of shape factors from c = 5.4 to c = 12.0. 
While c = 6.0 is the shape factor of a cube, each shape 
factor can correspond to a number of different grain shapes. 
The shape factor can be thought of as a deviation from the 
spherical case, the bigger the shape factor, the larger the 
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Figure 10. Nucleation rates as a function of time for six shape 
factors with 7 = 1.0. Filled circles indicate maximum nucleation 
rate. Rates shown are for nucleation within a shell of enclosed 
mass of ~ 4.96 to ~ 4.97 M Q . 



Figure 15. Total mass of carbon dust formed for all simulations. 



deviation from a sphere and the larger the surface for a 
given volume. 

In Figure [10] we show the nucleation rates for the same 
enclosed mass shell as shown in Figure [4] for increasing 
shape factors with 7 = 1.0. Since the shape factor appears 
in the exponential term of the nucleation rate equation as 
J s oc exp (— c 3 ) (Eq. (O), the higher shape factors reduce 
the nucleation rate. Therefore, the higher supersaturation 
levels at later times need to be attained so that the nucle- 
ation rate becomes high enough for significant depletion of 
the gas due to the growth of newly formed grains. However, 
the maximum nucleation rates for increasing shape factors 
are decreased. This is due to the fact that with a larger sur- 
face to volume ratio aspherical grains tend to grow faster, 
producing a sizeable depletion of the gas-phase carbon even 
for moderate nucleation rates. For all mass coordinates, in- 
creasing the shape factor leads to lower nucleation rates and 
delayed condensation times (Figures 1 1 lb and 112b . respec- 
tively). 

In Figure lllb -d we show the nucleation rate maxima 
for all shape factors for each of the other three sticking co- 
efficients. For any sticking coefficient, the shape factor has 
the same effect of suppressing the nucleation rate. All con- 
densation times, shown in Figure 112b — d. increase with the 
increased shape factor. Furthermore, for reduced sticking co- 
efficients, greater shape factors lead to greater delays in the 
condensation time. 

The shape factor c = (367t) 1,/3 is the only shape factor 
that has just one associated shape, the sphere. The other 
shape factors we consider do not necessarily have a unique 
grain shape. To find a size distribution that is easily com- 
parable to the spherical case, we define a volume equivalent 
radius for the aspherical grains as: 

/3V;\ 1/3 



V 4V J 



(6) 



where V is the grain volume. 

Figures [T3k -d show the size distributions for all shape 
factors and sticking coefficients. As the shape factor is in- 
creased, the size and number of larger grains also increases. 



Table 2. Total mass of carbon grains formed 



Mass (10" 2 M ) 



c 


7 = 1.0 


7 = 0.1 


7 = 0.01 


7 = 0.001 


(36tt) 1 /3 


4.95 


4.93 


4.89 


4.79 


5.4 


4.94 


4.91 


4.87 


4.75 


6.0 


4.92 


4.89 


4.84 


4.70 


7.0 


4.89 


4.86 


4.79 


4.62 


9.0 


4.82 


4.76 


4.65 


4.34 


12.0 


4.64 


4.52 


4.27 


3.32 



We find that the largest maximum grain radii are formed 
with c = 12.0 and 7 = 1.0, and have a volume equivalent ra- 
dius of almost 18 /im. On the other hand, the smallest maxi- 
mum grain radii (~ 0.008 /im) are formed with c = (367I") 1 / 3 
and 7 = 0.001. As the sticking coefficient is reduced, the 
maximum grain radii for each shape factor are also reduced, 
and the number of smaller grains is increased. Even though 
there is a much larger number of small grains than large 
grains, the majority of the dust mass is contained within 
intermediate sized grains. In Figure \T4h (7 = 1.0), most of 
the dust mass is contained in grains with volume equivalent 
radii between 0.01 and 0.5 /im. 

The masses of carbon grains formed (Figures 1 14b — d^) 
are dominated by the relatively small numbers of larger sized 
grains as the shape factor is increased. However, as the stick- 
ing coefficient is reduced (Figures I14b -d). even the masses 
of the grains formed with the largest shape factor become 
dominated by the smaller sized grains. The total mass of 
carbon grains formed for increased shape factors is given in 
Tableland shown in Figure IT51 Again, we notice that de- 
spite the difference in the size distribution due to different 
microphysical parameters, the total mass of dust that con- 
denses is fairly robust (within a factor of 1.5) and constitutes 
almost the total amount of carbon that was left in the gas 
phase after the creation of CO molecules. 
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Figure 11. Maximum nucleation rates for all six shape factors and four sticking coefficients. 
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Figure 17. Ry values for 23 simulations. Not shown is Ry = 41 
for c = 12.0 and 7 = 1.0. 



5 DISCUSSION 

We have shown that varying the microphysical properties 
of dust grains has important effects on the condensation 
times, nucleation rates, and size distributions of carbon dust 
grains from type II supernova explosions. However, the total 
mass of dust is only modestly affected by the changes in 
the grain properties. An inadequate choice of the shape or 
sticking coefficient is not therefore a viable explanation for 



the discrepancy between the mass of dust grains predicted 
in SN explosions and the observed dust mass in local type 
II SNe (see Section 1 for a more thorough discussion and 
references). 

We find that a larger saturation is necessary to achieve 
efficient nucleation with either a small sticking coefficient 
(7 < 1) or for aspherical grains (c > (36tt) 1/3 ). For that rea- 
son, all our simulations show that the condensation time 
grows when sticking coefficients less than unity or grain 
shapes that are different from spherical are adopted. How- 
ever, differences can be found in the nucleation rates and 
final size distribution of the grains. Simulations with a low 
value of the sticking coefficient show a delayed nucleation 
but very high nucleation rates, thereby producing large 
quantities of small grains. This is due to the fact that low 
sticking coefficients inhibit both nucleation and grain growth 
and, therefore, all the carbon remains in the gas phase un- 
til a high saturation level is reached. At that point, many 
small grai ns are nucleated and t he at omic carbon is quickly 
depleted. iBianchi fc Schneiderl (|2007l ) found similar effects 
when calculating dust nucleation with a sticking coefficient 
of 7 = 0.1. Asphericity of the grains, on the other hand, 
inhibits nucleation but enhances grain growth. As a conse- 
quence, even if fewer grains are nucleated, they grow fast 
and the result is a grain size distribution characterized by 
less numerous, larger grains. 

We find that the total mass of carbonaceous dust 
formed remains relatively stable even with sticking coeffi- 
cients as low as 0.001. For the spherical case only, we ex- 
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Figure 12. Condensation times for six shape factors and four sticking coefficients. 
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Figure 18. Total mass of dust formed for spherical carbonaceous 
grains with sticking coefficients down to 7 = 10 -9 . 



plored the possibility of even smaller sticking coefficients, 
down to 7 = 10 -9 . We find that the mass of carbon dust 
formed becomes significantly reduced for sticking coefficients 
of 7 = 10 -7 and below (see Figure fig)) . With sufficiently low 
values of sticking coefficient (below 10 -8 ) there is virtually 
no dust formation, but the required sticking coefficients seem 
unphysically low. 

In terms of the observable properties of the SN- 
condensed dust, we find that the quantity that is most af- 



fected is the extinction curve (see Figure I16J1 . Not surpris- 
ingly, simulations with small sticking coefficient (which, as 
explained above, produce large amounts of small grains), re- 
sult in a very steep extinction curve at far-UV wavelengths, 
with Kv values between 3 and 3.5, shown in Figure [T71 We 
make special note that the Kv values for 7 = 1.0 are rela- 
tively high, even though the size distribution of the grains 
is consistent with that of interstellar grains (see Figure |SJ) , 
because we account for only carbon grains and that includ- 
ing other grains, such as silicates, could decrease the Kv 
values. On the other hand, simulations with very aspherical 
grains and relatively high sticking coefficients produce larger 
grains and, consequently, grey extinction curves. Without 
the knowledge of the values of the sticking coefficient and of 
the shape factor it is therefore impossible to predict the ex- 
tinction curve of SN-condensed dust. This is a particularly 
worrying conclusion since the extinction curve is relatively 
easy to measure, even at high redshift, and could be used as 
an observational constraint for t he origin of dus t in th e vari- 
ous environments. For example, Maiol ino et alj ([2004) com- 
pared the extinction curve measured in a quasar at z = 6.2 
to the extinction cu r ve ca lculated using the dust model 
by iTodini fc Ferraral (|200lh . They find that the data and 
the theoretical prediction are in good agreement and con- 
clude that the dust observed in SDSSJ104845 05+463718.3 
is inde ed condensed in SN explosions (see also IStratta et al.l 
|2007h ). In light of our results, such conclusions need con- 
firmation once a complete theory of SN dust nucleation is 
obtained. 
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Figure 13. Size distribution of carbon grains for six shape factors and four sticking coefficients. 



The dust that condenses, however, is not the dust that 
is ejected into the ISM. Dust produced in a CCSN has to 
travel through the reverse shock before being released into 
the ISM. The reverse shock can destroy most of the dust, 
in pa r ticular the smaller dust grains (|Dra inc & Sa lpeterl 
19791: iNozawa et all [20071. iNath. Laskar. fc Shulll l200Sl : 
20101 ) . If astrophysical dust forma- 



Silvia. Smith, fc Shulll 



tion is indeed characterized by small values of the stick- 
ing coefficient, it is likely that the amount of dust formed 
is reduced significantly by the reverse shock. Conversely, in- 
creased shape factors allow for the formation of larger grains 
which would survive shock processing. The microphysical 
properties of dust grains can therefore affect the mass of dust 
that is injected in the interstellar medium, even though they 
affect only marginally the dust that is condensed during the 
early stages of the explosion. 

Grain nucleation with non-spherical shapes may be 
more complicated than we considered here. We have as- 
sumed that the shape factors of the grains do not change 
as the grains grow. Since we assume here that grains grow 
by the addition of monomers, the shape of small clusters 
can change as monomers attach, in turn altering the shape 
factor. Another route to take could be to nucleate grains at 
an arbitrary shape factor and then allow the grains to grow 
into spheres. For example, one may turn grains into spheres 
when the number of monomers is larger than a given value. 
In this case, increased shape factors may not lead to such 
large grains as we show in our size distributions (Figure PHI)) . 

This work is purely a parametric study that aims to 



show that the microphysical properties of the grains are im- 
portant, but not to point to any specific values of 7 and 
c to be used in nucleation calculations. Therefore, we have 
chosen to neglect important factors that need be considered 
in a complete SN dust nucleation model. These include the 
choi ce of progenitor model , the presence of other dust species 
(see |N ozawa et all (|2003h ) as well as charged molecules 
that may interfere with carbonaceous dust condensation, 
and the destruction of CO molecules, due to photodisso- 
ciation or collisons with fast electrons and charged particles 
ilPetuchowski et al.lll989l; iLepp. Dalgarno. fc McCravlll990l; 
Liu. Dalgarno. fc Lepd Il992l ; IClavton. Deneault. fc Mevej 
200 j) . that can inject additional c arbon atoms into the avai l- 



able monomer concentratio n (see iTodini fc Ferraral (|200ll ) ; 
iBianchi fc Schneiderl (|2007T )1. 

A zero-metallicity 20 Mq CCSN progenitor model was 
chosen because SNe at high redshift are expected to have 
zero metallicity. In general, the relative abundances of major 
elements in the He core are not significantly different among 
the SN models with different metallicities. Thus, the species 
of dust formed do not depend on th e metallicity of the SN 
progenitor star (|Nozawa et al.|[2010h . Additionally, the gas 
density and temperature in the He core are almost indepen- 
dent of the progenitor mass and metallicity as long as the 
kinet ic energy of the explosion is the same (N ozawa et alj 
2003). Therefore, SN models with non-zero metallicities, or 
with different progenitor masses, are expected to show sim- 
ilar effects on carbon grains formation as we see here. This 
means that the mass of carbon dust formed in the SN ejecta 
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Figure 14. Mass distribution of carbon grain for six shape factors and four sticking coefficients. 



is rather insensitive to changes in the sticking coefficient and 
shape factor and is purely determined by the mass of carbon 
atoms available for dust formation in the He layer. However, 
to confirm such expectations, additional progenitor models 
would need to be investigated. 

Perhaps more important than the choice in progenitor 
model is the dissociation process of CO molecules. In 
this paper we assumed the formation of CO molecules 
to be complete and considered only the condensa- 
tion process of C grains in the He layer where C/O 
> 1. In the expanding ejecta, CO molecules could 
be destroyed through interactions with fast electrons 
from radioactively decaying 56 Co and charged parti- 



cles such as He + and Ne + ([Petuchowski et al. 1989; 
Lepp. Dalgarno. fc McCrayl 199fj|: iLiu. Dalgarno. fc Leppl 
19921; IClavton. Liu, fc Dalgarnol 1 19991; 



Kwong, Chen, fc Fand l200d: IClayton. Deneault. fc MeveJ 
20011 ; iDeneault. Clayton, fc MeverlboOrj ). allowing for more 
free carbon (and oxygen) atoms to be available for grain 
formation than we consider here. However, the number 
abundance of silicon atoms is too small for most of the 
enclosed mass regions (M= 4.93-6.21 Mq) in this work 
(see Figure [3}, so that the formation of SiC and silicate 
gains cannot be expected. Therefore, the dissociation of 
CO molecules due to interactions with He + only results 
in a slight enhancement of the final mass of carbon grains 
and never changes our conclusion on the dependence 
of formation process of C grains on the microphysical 
properties. 



On the other hand, Clayton. Liu, fc Dalgarnol (|l999h 
and IDeneault. Clayton, fc Meverl i|2006T ) show that CO dis- 
sociation enables carbon dust grains to form even in O-rich 
layer where C/O < 1. Given that the abundance of silicon 
atoms in the O-rich layer is higher than in the He layer, the 
formation of S iC grains could be e xpected there. However, 
as discussed in lNozawa et all (12003), even if free carbon and 
silicon atoms coexist abundantly, the nucleation theory does 
not predict the formation of SiC grains. The formation pro- 
cess of large SiC grains as appeared in presolar grains, as 
well as formation efficiency of molecules is to be pursued in 
more sophisticated studies of dust formation. Furthermore, 
in the O-rich layer, the formation o f silicate grains is also 
feasible. iBianchi fc Schneiderl i|2007l ) show that the forma- 
tion of silicate grains is more sensitive to changes in sticking 
coefficient than carbon grains, and the mass of silicate grains 
formed can be reduced for even 7 = 0.1. The inclusion of 
silicate grains could affect the resulting extinction curves. 



We have adopted the thermodynamic approach for this 
study because it involves the simplest nucleation equations, 
however, use of the kinetic theory of nucleation should be 
considered in the future. The kinetic theory still needs to 
take the sticking coefficient into account, but makes consid- 
eration of an evolving shape factor unnecessary, because the 
shape of the grain from a complex solid (for a few molecules) 
to a sphere (for ~ 100 molecules) is intrinsically taken into 
account. 
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Figure 16. Extinction curves for each sticking coefficient and shape factor. 
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